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Severe acute respiratory syndrome (SARS) emerged in 2002 to 2003 in southern 
China. The origin of its etiological agent, the SARS coronavirus (SARS-CoV), 
remains elusive. Here we report that species of bats are a natural host of 
coronaviruses closely related to those responsible for the SARS outbreak. These 
viruses, termed SARS-like coronaviruses (SL-CoVs), display greater genetic 
variation than SARS-CoV isolated from humans or from civets. The human and 
civet isolates of SARS-CoV nestle phylogenetically within the spectrum of SL- 
CoVs, indicating that the virus responsible for the SARS outbreak was a member 


of this coronavirus group. 


Severe acute respiratory syndrome (SARS) 
was caused by a newly emerged coronavirus, 
now known as SARS coronavirus (SARS-CoV) 
(, 2). In spite of the early success of etio- 
logical studies and molecular characterization 
of this virus (3, 4), efforts to identify the ori- 
gin of SARS-CoV have been less successful. 
Without knowledge of the reservoir host dis- 
tribution and transmission routes of SARS- 
CoV, it will be difficult to prevent and control 
future outbreaks of SARS. 

Studies conducted previously on animals 
sampled from live animal markets in Guang- 
dong, China, indicated that masked palm civets 
(Paguma larvata) and two other species had 
been infected by SARS-CoV (5). This led to 
a large-scale culling of civets to prevent fur- 
ther SARS outbreaks. However, subsequent 
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studies have revealed no widespread infection 
in wild or farmed civets (6, 7). Experimental 
infection of civets with two different human 
isolates of SARS-CoV resulted in overt clin- 
ical symptoms, rendering them unlikely to be 
the natural reservoir hosts (8). These data sug- 
gest that although P. larvata may have been 
the source of the human infection that pre- 
cipitated the SARS outbreak, infection in this 
and other common species in animal markets 
was more likely a reflection of an “artificial” 
market cycle in naive species than an indi- 
cation of the natural reservoir of the virus. 

Bats are reservoir hosts of several zoonotic 
viruses, including the Hendra and Nipah viruses, 
which have recently emerged in Australia and 
East Asia, respectively (9-11). Bats may be 
persistently infected with many viruses but rarely 
display clinical symptoms (/2). These charac- 
teristics and the increasing presence of bats and 
bat products in food and traditional medicine 
markets in southern China and elsewhere in 
Asia (/3) led us to survey bats in the search 
for the natural reservoir of SARS-CoV. 

In this study, conducted from March to 
December of 2004, we sampled 408 bats 
representing nine species, six genera, and 
three families from four locations in China 
(Guangdong, Guangxi, Hubei, and Tianjin) af- 
ter trapping them in their native habitat (Table 
1). Blood, fecal, and throat swabs were col- 
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lected; serum samples and cDNA from fecal 
or throat samples were independently analyzed, 
double-blind, with different methods in our 
laboratories in Wuhan and Geelong (/4). 

Among six genera of bat species surveyed 
(Rousettus, Cynopterus, Myotis, Rhinolophus, 
Nyctalus, and Miniopterus), three communal, 
cave-dwelling species from the genus Rhinolo- 
phus (horseshoe bats) in the family Rhino- 
lophidae demonstrated a high SARS-CoV 
antibody prevalence: 13 out of 46 bats (28%) 
in R. pearsoni from Guangxi, 2 out of 6 bats 
(33%) in R. pussilus from Guangxi, and 5 out 
of 7 bats (71%) in R. macrotis from Hubei. 
The high seroprevalence and wide distribution 
of seropositive bats is expected for a wildlife 
reservoir host for a pathogen (/5). 

The serological findings were corroborated 
by poylmerase chain reaction (PCR) analyses 
with primer pairs derived from the nucleocapsid 
(N) and polymerase (P) genes (table S1). Five 
fecal samples tested positive, all of them from 
the genus Rhinolophus: three in R. pearsoni 
from Guangxi and one each in R. macrotis and 
R. ferrumequinum, respectively, from Hubei. 
No virus was isolated from an inoculation of 
Vero E6 cells with fecal swabs of PCR-positive 
samples. 

A complete genome sequence was deter- 
mined directly from PCR products from one of 
the fecal samples (sample Rp3) that contained 
relatively high levels of genetic material. The 
genome organization of this virus (Fig. 1), ten- 
tatively named SARS-like coronavirus isolate 
Rp3 (SL-CoV Rp3), was essentially identical to 
that of SARS-CoV, with the exception of three 
regions (Fig. 1, shaded boxes). The overall 
nucleotide sequence identity between SL-CoV 
Rp3 and SARS-CoV Tor2 was 92% and in- 
creased to ~94% when the three variable re- 
gions were excluded. The variable regions are 
located at the 5’ end of the S gene (equivalent to 
the S1 coding region of coronavirus S protein) 
and the region immediately upstream of the 
N gene. These regions have been identified 
as “high mutation” regions among different 
SARS-CoVs (5, 16, 17). The region upstream 
of the N gene is known to be prone to de- 
letions of various sizes (5, 16, 18). 

Predicted protein products from each gene 
or putative open reading frame (ORF) of SL- 
CoV Rp3 and SARS-CoV Tor2 were com- 
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Table 1. Detection of antibodies to SARS-CoV and PCR amplification of N and P gene fragments with 
SARS-CoV-specific primers. ND, not determined because of poor sample quality or unavailability of 


specimens from individual animals. 


? PCR analysis: 
Sampling . Antibody test: positive/total (%) 
Babe specles positive/total (%) 

. P Fecal Respiratory 
Time Location ewabs swabs 
Mar 04 —— Nanning, Guangxi Rousettus leschenaulti 1/84 (1.2%) 0/110 ND 

Maoming, Guangdong _—Rousettus leschenaulti 0/42 0/45 ND 
Cynopterus sphinx 0/17 0/27 ND 

July 04 ~~ Nanning, Guangxi Rousettus leschenaulti ND 0/55 0/55 
Tianjin Myotis ricketti ND 0/21 0/21 

Nov 04 Yichang, Hubei Rhinolophus pusillus ND 0/15 ND 
Rhinolophus ferrumequinum 0/4 1/8 (12.5%) * ND 

Rhinolophus macrotis 5/7 (71%) —- 1/8 (12.5%) + 0/3 

Nyctalus plancyi 0/1 0/1 ND 

Miniopterus schreibersi 0/1 0/1 ND 

Myotis altarium 0/1 0/1 ND 

Dec 04 —— Nanning, Guangxi Rousettus leschenaulti 1/58 (1.8) ND ND 
Rhinolophus pearsoni 13/46 (28.3%) 3/30 (10%)t 0/11 

Rhinolophus pussilus 2/6 (33.3%) 0/6 0/2 


*Positive fecal sample designated Rf1 
designated Rp1, Rp2, and Rp3, respectively. 


A 


{Positive fecal sample designated Rm1 


{Positive fecal samples 
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Fig. 1. Genome organization of, and comparison between, SL-CoV and SARS-CoV. (A) Overall 
genome organization of SL-CoV Rp3. (B) Expanded diagram of the 3’ region of the genome in 
comparison with SARS-CoV strains Tor2 and SZ3, following the same nomenclature used by Marra 
et al. (4). The genes (named by letters P, S, E, M, and N) present in all coronaviruses are shown in 
dark-colored arrows, whereas the SARS-CoV-specific ORFs are numbered and illustrated in light- 
colored arrows. ORF10’ follows the nomenclature by Guan et al. (5) to indicate that the single ORF 
present between ORF9 and N in SL-CoV is equivalent to the fusion of ORF10 and ORF11 in the 
same region in SARS-CoV Tor2. The shaded boxes mark the only three regions displaying 
significant sequence difference between the two viruses (table $2). 


pared (table S2). The P, S, E, M, and N pro- 
teins, which are present in all coronaviruses, 
were similarly sized in the two viruses, with 
sequence identities ranging from 96% to 
100%. The only exception was the S1 domain 
of the S protein, where sequence identity fell to 
64%. The S1 domain is involved in receptor 
binding, whereas the S2 domain is responsible 
for the fusion of virus and host cell membranes 
(19). The sequence divergence in the S1 do- 
main corroborated our serum neutralization 
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studies, which indicated that although bat sera 
have a high level of cross-reactive antibodies 
(with enzyme-linked immunosorbent assay 
titers ranging from 1:100 to 1:6400), they 
failed to neutralize SARS-CoV when tested on 
Vero E6 cells. This finding suggests that S1 is 
the main target for antibody-mediated neutral- 
ization of this group of viruses, which is con- 
sistent with previous reports indicating that 
major SARS-CoV neutralization epitopes are 
located in the S1 region (20, 2/). 
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In addition to the five genes present in all 
coronavirus genomes, coronaviruses also have 
several ORFs between the P gene and the 3’ 
end of the genome that code for nonstructural 
proteins. The function of these nonstructural 
proteins is largely unknown. The location and 
sequence of ORFs are group- or virus-specific 
and hence can serve as important molecular 
markers for studying virus evolution and clas- 
sification (19, 22). SARS-CoV has a unique 
set of ORFs not shared by any of the known 
coronaviruses (3, 4). Most of these ORFs 
were also present in SL-CoV, confirming the 
extremely close genetic relationship between 
SARS-CoV and SL-CoV (Fig. | and table $2). 

Coronaviruses produce subgenomic mRNAs 
through a discontinuous transcription process 
not fully characterized (/9). Conserved nucle- 
otide sequences functioning as transcription 
regulatory sequences (TRSs) are required for 
the production of the subgenomic mRNAs. In 
SARS-CoV, such TRSs were identified at 
each of the predicted gene start sites (3, 4). 
All of these TRSs were absolutely conserved 
between SARS-CoV Tor2 and SL-CoV Rp3 
(table S3), further demonstrating that these 
two viruses are very closely related. 

SL-CoV is completely different from a 
bat coronavirus (bat-CoV) recently identified 
by Poon et al. (7) from species of bats in the 
genus Miniopterus during a wildlife surveil- 
lance study in Hong Kong (Fig. 2). Because 
the complete genome sequence was not 
available for bat-CoV, only the trees covering 
the common sequences (i.e., parts of the P1b 
and S2 proteins) are shown. The phylogenetic 
analysis demonstrated that SL-CoV Rp3 and 
SARS-CoVs are clustered together but that 
bat-CoV is placed among the relatively dis- 
tant group | viruses. Hereafter, SARS-CoVs 
and SL-CoVs will be collectively called the 
SARS cluster of coronaviruses. 

In addition to the complete genome se- 
quence of SL-CoV Rp3, partial genome se- 
quences for the other four PCR-positive bat 
samples were also determined. Phylogenetic 
analysis based on the N protein sequences 
(Fig. 3A) revealed that the genetic variation 
among the SL-CoV sequences was much greater 
than that exhibited by SARS-CoVs (for sim- 
plicity, only three human and civet SARS-CoV 
isolates were used; the remainder are almost 
identical to those shown). This was especially 
obvious when SL-CoVs isolated from different 
bat species were compared. Moreover, the re- 
sults suggested that SARS-CoVs nestle phylo- 
genetically within the spectrum of SL-CoVs. 

We also compared the “high mutation” re- 
gions in samples Rfl, Rm1, and Rp3. For the 
region upstream of the N gene, SL-CoVs from 
all three bat species contained a single ORF 
(ORF10’), similar to that found in SARS-CoV 
isolates from civets (5) and patients in the 
early phase of the outbreaks (/6, 18) but dif- 
ferent from that in most human isolates, which 
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Fig. 2. Phylogenetic $Z3 
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Fig. 3. Phylogenetic trees based on deduced ami- 
no acid sequences of (A) N, (B) ORF10’, and (C) 
S1 proteins. Tor2, $Z3, and GDO1, different SARS- 
CoV strains; Rf1, Rm1, and Rp1-3, different SL- 
CoV sequences. The genetic distance scale shown 
for (A) is different from those for (B) and (C). 


have a 29-nucleotide deletion in this region 
(3, 4, 16). ORF10’ in Rf1 codes for a protein 
having the same size (122 amino acids) as and 
more than 80% sequence identity to ORF10’ 
proteins of SARS-CoVs, but those in Rm1 
and Rp3 code for a 121—amino acid protein 
with only 35% sequence identity (Fig. 3B and 
fig. S2). By contrast, analysis of the S1 pro- 
tein regions (Fig. 3C and fig. S3) indicated 
that Rfl was more closely related to SL-CoVs 
from two other bat species than to SARS- 
CoVs, suggesting that the SARS cluster of 
coronaviruses could recombine to increase 
genetic diversity and fitness, as is well doc- 
umented for other coronaviruses (/9). We 
were unable to sequence these regions for Rp1 
or Rp2, owing to the poor quality of the fecal 


materials from these two animals. The limited 
amount of cDNA available was used up for N 
gene analysis and in initial sequencing trials 
with SARS-CoV-derived primers, which 
were largely unsuccessful. Judging from the 
close relationship of the N genes between 
Rp1, Rp2, and Rp3 (fig. $1), it is unlikely that 
Rp! or Rp2 will have major sequence differ- 
ences from Rp3 in the S1 or ORF10’ regions. 
This is not unexpected, considering that these 
three positive samples were obtained from the 
same bat species in the same location. 

The genetic diversity of bat-derived se- 
quences supports the notion that bats are a 
natural reservoir host of the SARS cluster of 
coronaviruses. A similar observation has been 
made for henipaviruses, another important group 
of emerging zoonotic viruses of bat origin, 
which show greater genetic diversity in bats 
than was observed among viruses isolated 
during the initial Nipah outbreaks in Malaysia 
(23-26). The overall nucleotide sequence iden- 
tity of 92% between SL-CoVs and SARS- 
CoVs is very similar to that observed between 
Nipah viruses isolated from Malaysia and 
Bangladesh in 1999 and 2004, respectively 
(25) (fig. S4). SL-CoVs present a new chal- 
lenge to the diagnosis and treatment of future 
disease outbreaks. The current tests and ther- 
apeutic strategies may not work effectively 
against all viruses in this group, owing to 
their great genetic variability in the S1 do- 
main region of the S gene. 

The genus Rhinolophus contains 69 species 
and has a wide distribution from Australia to 
Europe (27). They roost primarily in caves and 
feed mainly on moths and beetles. However, 
notwithstanding the predominant Rhinolophus 
findings in this study, it is highly likely that 
there are more SARS-related coronaviruses to 
be discovered in bats. Indeed, our positive 
serological findings in the cave-dwelling fruit 
bat Rousettus leschenaulti indicate that infec- 


tion by a related virus could occur in fruit bats 
as well, albeit at a much lower frequency. A 
plausible mechanism for emergence from a 
natural bat reservoir can be readily envisaged. 
Fruit bats including R. leschenaulti, and less 
frequently insectivorous bats, are found in mar- 
kets in southern China. An infectious consign- 
ment of bats serendipitously juxtaposed with a 
susceptible amplifying species, such as P. 
larvata, at some point in the wildlife supply 
chain could result in spillover and establish- 
ment of a market cycle while susceptible ani- 
mals are available to maintain infection. Further 
studies in field epidemiology, laboratory infec- 
tion, and receptor distribution and usage are 
being conducted to assess potential roles played 
by different bat species in SARS emergence. 
These findings on coronaviruses, together 
with data on henipaviruses (23-25, 28), sug- 
gest that genetic diversity exists among zoo- 
notic viruses in bats, increasing the possibility 
of variants crossing the species barrier and 
causing outbreaks of disease in human popula- 
tions. It is therefore essential that we en- 
hance our knowledge and understanding of 
reservoir host distribution, animal-animal and 
human-animal interaction (particularly within 
the wet-market system), and the genetic diversity 
of bat-borne viruses to prevent future outbreaks. 
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Ciliary neurotrophic factor (CNTF) induces weight loss in obese rodents and 
humans, and for reasons that are not understood, its effects persist after the 
cessation of treatment. Here we demonstrate that centrally administered CNTF 
induces cell proliferation in feeding centers of the murine hypothalamus. Many 
of the newborn cells express neuronal markers and show functional phenotypes 
relevant for energy-balance control, including a capacity for leptin-induced phos- 
phorylation of signal transducer and activator of transcription 3 (STAT3). Co- 
administration of the mitotic blocker cytosine-B-D-arabinofuranoside (Ara-C) 
eliminates the proliferation of neural cells and abrogates the long-term, but not 
the short-term, effect of CNTF on body weight. These findings link the sustained 
effect of CNTF on energy balance to hypothalamic neurogenesis and suggest that 
regulated hypothalamic neurogenesis in adult mice may play a previously 
unappreciated role in physiology and disease. 


The obesity epidemic has prompted major ef- 
forts to develop safe and effective therapies 
(1, 2). However, approved drugs for obesity 
have limited efficacy and act only acutely, 
with patients rapidly regaining weight after 
terminating treatment (3). Only the neuro- 
cytokine ciliary neurotrophic factor (CNTF) 
and Axokine, an analog of CNTF developed 
as a drug candidate for the treatment of obe- 
sity, appear to deviate from this paradigm. 
Rodents and patients treated with Axokine 
were reported to maintain lowered body 
weights weeks to months after the cessation 
of treatment (4, 5). This feature of Axokine/ 
CNTF action is unexplained and suggests that 
CNTF induces long-lasting changes in one or 
more elements of the energy-balance circuitry. 

In rodents, CNTF is most potent when 
administered directly into the cerebrospinal 
fluid (6) and activates signaling cascades in 
hypothalamic nuclei involved in feeding 
control (5, 7, 8). For instance, CNTF acti- 
vates phosphorylation of signal transducer 
and activator of transcription 3 (STAT3) in a 
population of hypothalamic neurons that 
substantially overlaps with those activated 
by leptin (5). However, in contrast to CNTF, 
leptin-treated animals do not maintain their 
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lowered body weight after the cessation of 
treatment. We thus sought a CNTF-specific 
mechanism to explain this long-term effect. 
CNTF supports the survival of neurons in 
vitro and in vivo (9) and has also been impli- 
cated in the maintenance of adult neural stem 
cells (J0). Furthermore, other trophic factors, 
such as epidermal growth factor and fibroblast 
growth factor 2, are known to act as mitogens 
on adult neuronal progenitors (1/7, /2), and they 
promote the functional regeneration of hippo- 
campal pyramidal neurons (/3). Neurogenesis 
in the adult brain is most clearly defined in the 
subventricular zone (SVZ) of the lateral ven- 
tricles and the subgranular zone (SGZ) of the 
hippocampal formation (/4). However, recent 
reports indicate that the neuroproliferative 
potency in the adult extends to other brain 
structures, including the hypothalamus (/5—/7). 
On the basis of these findings, we hypothesized 
that the long-term effect of CNTF on body- 
weight regulation might involve neurogenesis 
in the hypothalamus, which is the brain region 
most relevant for energy-balance regulation. 
To assess the mitogenic potency of CNTF 
in the adult nervous system in vivo, we 
delivered the cell-proliferation marker bromo- 
deoxyuridine (BrdU) alone (vehicle treatment) 
or together with CNTF directly into the 
cerebrospinal fluid of mouse brains (/8). 
CNTF and BrdU were continuously infused 
for 7 days into the right lateral ventricle using 
osmotic minipumps. Mice were switched to a 
high-fat diet two months before surgery and 
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were kept on this diet throughout the experi- 
ments. In accordance with previous results (5), 
CNTF-treated mice showed a marked reduc- 
tion in body weights (Fig. 1A), which persisted 
after termination of CNTF delivery. Mice were 
killed 22 days after surgery, and brain sections 
were immunostained with an antibody against 
BrdU. Because BrdU incorporates into DNA of 
dividing cells, BrdU-positive (BrdU*) cells are 
thought to represent newborn cells. Figure 1B 
shows coronal sections of vehicle- and CNTF- 
infused animals at the level of the arcuate, 
ventromedial, and dorsomedial nuclei, well- 
known hypothalamic centers for energy-balance 
regulation (/9). In vehicle-infused animals, few 
BrdUt cells were detected in the parenchyma 
surrounding the third ventricle (Fig. 1B, left). 
Administration with CNTF led to a dramatic 
increase of BrdU+ cells (Fig. 1B, right). Note 
the higher density of BrdU* cells at the base 
of the third ventricle, which is part of the ar- 
cuate nucleus/median eminence. 

The pattern of CNTF receptor (CNTFR) 
mRNA expression is consistent with this ob- 
servation. In situ hybridization using a ribo- 
probe against CNTFR mRNA revealed strong 
staining in the walls of the basal third ventricle 
and surrounding arcuate nucleus parenchyma 
(Fig. 1C). Because this section originated from 
an animal treated with both CNTF and BrdU, 
we colabeled with antibodies to BrdU. Many 
BrdU* cells were positive for CNTFR ex- 
pression, indicating that CNTF, at least in part, 
directly promotes cell division by binding to 
CNTFR on putative neural progenitor cells 
(Fig. 1D, inset). By counting all newly gener- 
ated cells in the caudal hypothalamus, CNTF 
treatment led to a marked increase of BrdUt 
cells over vehicle-infused animals (Fig. 1E). 
The total number of BrdU* cells in CNTF- 
treated animals remained constant for at least 
2 weeks after the infusion period. Subsequent- 
ly, the numbers decreased but plateaued at a 
high level. Vehicle-infused animals showed a 
similar fractional decrease over time. Thus it 
appears that the majority of hypothalamic 
BrdU* cells do not die or migrate to distant 
areas as reported for newborn neurons of the 
SVZ, which follow the rostral migratory stream 
toward the olfactory bulb (20). 

To investigate the origin of adult-born cells 
in the hypothalamus, we examined CNTF and 
vehicle-infused brains every 12 hours starting 
48 hours after surgery, a time when the in- 
fused CNTF/BrdU should just reach the ven- 
tricular system (/8). Hypothalamic BrdU 
incorporation was first detected 60 to 72 hours 
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